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Abstract 

We present results of dynamical simulations of Nf = 2 degenerate Wilson 
twisted mass quarks at maximal twist in the range of pseudo scalar masses 
300 MeV < mps < 550 MeV. Reaching such small masses was made pos- 
sible owing to a recently developed variant of the HMC algorithm. The 
simulations are performed at one value of the lattice spacing a < 0.1 fm. In 
order to have 0(a) improvement and aiming at small residual 0(a?) cutoff 
effects, the theory is tuned to maximal twist by requiring the vanishing of 
the untwisted quark mass. Precise results for the pseudo scalar decay con- 
stant and the pseudo scalar mass are confronted with chiral perturbation 
theory predictions and the low energy constants F, I3 and Z4 are evaluated 
with small statistical errors. 



1 Introduction 

The Wilson twisted mass formulation of lattice QCD, though a rather recent 
approach, has been by now well established. It amounts to adding a twisted 
mass term to the standard, un-improved Wilson-Dirac operator [1] leading to 
so-called Wilson twisted mass fermions [21 [3] . 

Besides being a theoretically sound formulation of lattice QCD, Wilson 
twisted mass fermions offer a number of advantages when tuned to maximal 
twist: (i) in this case automatic 0(a) improvement is obtained by tuning only 
one parameter, the bare untwisted quark mass, while avoiding additional tuning 
of operator-specific improvement-coefficients; (ii) the mixing pattern in the renor- 
malisation process can be significantly simplified; (iii) the twisted mass provides 
an infra-red regulator helping to overcome possible problems with ergodicity in 
molecular dynamics based algorithms^. 

In the quenched approximation, these expectations - based on general field 
theoretical and chiral perturbation theory (xPT) related arguments [21 El HH [121 
[131 IH US] - could be verified in actual simulations PSI QU QSl ESJ: 0(a) im- 
provement is indeed realised when the theory is tuned to maximal twist. More- 
over, it has been shown that a particular realisation of maximal twist, requiring 
parity restoration, also suppresses the 0(a 2 ) cut-off effects substantially, even 
at small quark masses corresponding to values of the pseudo scalar mass of 
mps < 300MeV. In addition, with the twisted mass parameter as an infra-red 
cut-off in place, substantially smaller quark masses could be obtained, compared 
to those reachable by standard or 0(a) improved Wilson fermions which are 
plagued by so-called exceptional configuration problems in the quenched approx- 
imation. In Refs. [2ni EU E2] it was shown that "wrong chirality" mixing effects 

1 Although in the light of recent algorithmic developments [31 [5l El HI [9] this property 
does not seem to be that important anymore, we consider it still to be an advantage to have 
an infra-red regulator in the theory which helps in stabilising the simulations. For a recent 
stability analysis of pure Wilson fermion simulations see Ref. [TU] 
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in the renormalisation process are substantially reduced. In Ref. [TT] it was 
proved that all such mixings can be eliminated if a mixed action with maximally 
twisted sea quarks and appropriately chosen Osterwalder-Seiler valence fermions 
is employed. For a further discussion of the potential of twisted mass QCD on 
the lattice, see Refs. [28| [2^ [25| [26] . 

The main drawback of the twisted mass approach is the explicit breaking of 
parity and isospin symmetry which are only restored when the continuum limit is 
reached. However, due to 0(a) improvement, this breaking is an 0(a 2 ) effect as 
confirmed by simulations performed in the quenched approximation [271 EE] • For 
recent reviews of the status of Wilson twisted mass fermions see Refs. [291 GO, ED] 
and references therein. 

It is the main goal of our collaboration to compute a number of phenomeno- 
logically relevant quantities with dynamical quarks, in the continuum limit and 
at small values of the pseudo scalar mass. As a first step in this direction we 
here present results for Nf = 2 mass-degenerate quarks at a fixed lattice spacing 
a < 0.1 fm. We have so far concentrated on the pseudo scalar mass mpg, covering 
a range of values 300MeV < mpg < 550MeV, the pseudo scalar decay constant 
fps and the static inter-quark force parameter r at five values of the quark mass. 
A wider range of physical observables will be addressed in the future. The results 
for mps and /pg are confronted with predictions of %PT which allows extracting 
the low-energy constants Z3, U, F and Bq of the corresponding effective chiral 
Lagrangian. We also provide a determination of the size of isospin violation mea- 
sured from the mass splitting between the lightest charged and neutral pseudo 
scalar particles. First accounts of our work were presented at recent conferences, 
see Refs. j32j [33] . In this publication we will focus on the results of our present 
simulations obtained at one value of (3 and one volume. We shall provide, in 
a forthcoming paper [M], a comprehensive description of our analysis procedure 
and address systematic errors by including future runs on larger lattices, at differ- 
ent values of (3 and with extended statistics. Related works with Wilson fermions 
at similar small pseudo scalar meson masses are published in Refs. [351 1361 [37] . 

2 Choice of Lattice Action 

The Wilson twisted mass fermionic lattice action for two flavours of degenerate 
quarks reads (in the so called twisted basis [2] and fermion fields with continuum 
dimensions) 
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where am is the bare untwisted quark mass and an the bare twisted mass, T3 
is the third Pauli matrix acting in flavour space and r is the Wilson parameter, 
which we set to one in our simulations. Twisted mass fermions are said to be at 
maximal twist if the bare untwisted mass is tuned to its critical value, m cr it- We 
will discuss later how this can be achieved in practice. 

In the gauge sector we use the so called tree-level Symanzik improved gauge 
action (tlSym) [3B], which includes besides the plaquette term also rectan- 

gular (1x2) Wilson loops U^ 2 U 

s s = f E ( 6 ° E - ReTr (C)}+ 6 i E - ReTr (^)} ) ( 2 ) 

X \ /i,i/ = l /i,i/ = l / 

with /? the bare inverse coupling, b\ = —1/12 and the (proper) normalisation 
condition bo = 1 — 8bi. Note that at b\ = this action becomes the usual Wilson 
plaquette gauge action. 



2.1 0(a) improvement 

As mentioned before, 0(a) improvement can be obtained by tuning Wilson 
twisted mass fermions to maximal twist. In fact, it was first proved in Ref. [3] that 
parity even correlators are free from 0(a) lattice artifacts at maximal twist by 
using spurionic symmetries of the lattice action. Later on it was realised (T2J [31] 
that a simpler proof is possible based on the parity symmetry of the continuum 
QCD action and the use of the Symanzik effective theory. 

From this latter way of proving 0(a) improvement, it becomes also clear 
how to define maximal twist: first, choose an operator odd under parity (in the 
physical basis) which has a zero expectation value in the continuum. Second, 
at a non-vanishing value of the lattice spacing tune the expectation value of 
this operator to zero by adjusting the value of am . This procedure, which has 
been proposed in [I2B ED] and has been theoretically analysed in [12] , is sufficient 
to define maximal twist independently of the chosen operator. To approach 
smoothly the continuum limit this tuning has to be performed at fixed physical 
situation while decreasing the lattice spacing. 

It was shown in an extended scaling test in the quenched approximation, 
that 0(a) improvement works extremely well for maximally twisted mass quarks 
[Ttjj [T71 HE]- In the context of this scaling test, the method of setting the so-called 
PCAC mass to zero was found to be very successful, in agreement with theoretical 
considerations [T5 l [TBI fT2] . Here the PCAC mass 

mpcAc- 2Ex(pfl(l)pa(0)) > a" 1 ' 2 ( 3 ) 

is evaluated at large enough time separation, such that the pion ground state 
is dominant. To see that the procedure of defining am cr - lt from the vanishing of 
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^pcac is the appropriate one, it is enough to recall that under that condition 
the multilocal operator ^ x AQ(x)P a (0) becomes, in the physical basis, the parity 
odd operator e 3ab J2xi J lo'r b ip(x) ip^T a ip{Q). 

In principle one could think of determining om cr i t at each value of a/i at which 
simulations are performed, possibly followed by an extrapolation to vanishing 
a/i. The strategy we are following in this paper is, instead, to take the value of 
am CT [ t from the simulation at the lowest available value a/x m ; n <C aAqcD- In this 
situation 0(a) improvement is still guaranteed, because working at /i min merely 
leads to 0(a/x m i n AQCD) effects in m crit and C(a 2 / u m i n AQCD) relative corrections in 
physical quantities [12] . 

2.2 Phase Structure 

In order to understand our choice of the gauge action, it is important to realise 
that Wilson-type fermions have a non-trivial phase structure at finite lattice spac- 
ing: in a series of publications [HI [391 SO, [291 S21 S3] the phase structure of lattice 
QCD was explored. For lattice spacings a > 0.15 fm clear signals of first order 
phase transitions at the chiral point were found. The strength of those phase 
transitions weakens when the continuum limit is approached. This phase tran- 
sition was identified to be a generic property of Wilson-type fermions since the 
phenomenon takes place for the pure Wilson as well as the Wilson twisted mass 
formulation [21 E] of lattice QCD. Also the properties of physical quantities mea- 
sured in both metastable branches of this first order phase transition were studied 
and compared to results of xPT [HI HEl EH SSJ ES IE] finding that (lattice) %PT 
describes the simulation data quite well. This is somewhat surprising since the 
simulation data were obtained at rather coarse values of the lattice spacing and 
at rather heavy pseudo scalar masses, where the applicability of ;\RT may be 
questionable. 

A very important consequence of the first order phase transition phenomenon 
is that at non-vanishing lattice spacing, simulations cannot be performed with 
pseudo scalar mesons below a certain minimal mass value. From lattice 
analyses it is expected that this minimal value of the pseudo scalar mass goes to 
zero with a rate of 0(a). In different words, given a value of the pseudo scalar 
mass, m PS , one can always find a value of the lattice spacing a max (mps), such 
that simulations at a < a max (mps) can be safely performed. For example, when 
the Wilson plaquette gauge action is used one finds a max ~ 0.07 fm to realise a 
pseudo scalar mass of about 300MeV [4"3] . 

The phase structure of lattice QCD with Wilson-type fermions has previously 
been addressed: there have been investigations concerning the Aoki-phase [IB] 
in Refs. [HU [50] at large gauge couplings corresponding to values of the lattice 
spacing a > 0.2 fm. In other studies [ED, [52], [53] signals of first order phase 
transitions were found for Wilson fermions with and without the clover term, see 
also Ref . [51] . In Ref . [53] a speculative picture of the phase structure of Wilson 
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L 3 • T 




^crit (c/^min) 


r /a 


3.9 


24 3 • 48 


0.004 


0.160856 


5.22(2) 



Table 1: Simulation parameters. We denote by a/x m ; n the smallest value of the twisted mass 
parameter a[i at which we have performed simulations. At this value of a[i we determined the 
critical mass m cr ;t, or, equivalently the critical hopping parameter K cr j t = 1/(8 + 2om cr it). The 
value of r /a has been extrapolated to the physical point, where tops = 139.6 MeV. 



lattice QCD has been given and in Ref. [56] an analysis within the framework of 
%PT has been reported. A detailed understanding of the generic phase structure 
was obtained in the 2-dimensional Gross-Neveu model, see Refs. [57J EEl EH] • Of 
course, it is unclear how much these last results are applicable to 4-dimensional 
lattice QCD. 

In order to choose a gauge action for our production simulations we studied 
the phase structure employing a number of different gauge actions: the stan- 
dard Wilson plaquette gauge action pQ (&i = in Eq. (j2J)), the DBW2 gauge 
action [BU] ipi = —1.4088) and the tree-level Symanzik improved gauge action 
[38] (pi = —1/12). A marked dependence of the strength of the phase transition 
on the choice of the gauge action has been found. In particular, these investiga- 
tions revealed that the DBW2 and the tlSym gauge actions substantially weaken 
the effect of the first order phase transition and in particular the value of a max 
increases when the coefficient b\ in Eq. is moved away from zero [101 H3] . We 
refer to Refs. [29l [31] for summaries of these results. 

The DBW2 gauge action appears to lead to a bad scaling behaviour [6T 1 1621 163] 
and a slow convergence of perturbation theory [[64], whereas the tlSym gauge 
action is expected to show - by construction - a good scaling behaviour and a 
fast convergence of perturbation theory. Therefore, the tlSym gauge action looks 
like a good compromise between the Wilson gauge action which is most strongly 
affected by the first order phase transition and the DBW2 gauge action. 

3 Numerical Results 
3.1 Set-up 

In this letter we will present results at a fixed value of the lattice spacing of 
a < 0.1 fm only. In tabled] we provide the value of a/x m i n at which we imposed 
the vanishing of mpcAc (Eq. (G2)) and thus determined om crit . In table [2] we list 
the values of the quark mass ampcACi the pseudo scalar mass amps, the pseudo 
scalar decay constant a/ps, r /a and the plaquette integrated autocorrelation 
time at all values of the twisted mass parameter a\i. All other parameters were 
kept fixed as specified in table [H 

The algorithm we used is a HMC algorithm with mass preconditioning [U [65] 
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a/i 


amps 


ajps 


ampcAC 


r /a 


Tmt{P> 


0.0040 


0.13587(68) 


0.06531(40) 


-0.00001(27) 


5.196(28) 


55(17) 


0.0064 


0.16937(36) 


0.07051(35) 


-0.00009(17) 


5.216(27) 


23(07) 


0.0085 


0.19403(50) 


0.07420(24) 


-0.00052(17) 


5.130(28) 


13(03) 


0.0100 


0.21004(52) 


0.07591(40) 


-0.00097(26) 


5.143(25) 


15(04) 


0.0150 


0.25864(70) 


0.08307(34) 


-0.00145(42) 


5.038(24) 


06(02) 



Table 2: Results from simulations at (3 = 3.9 using the simulation parameters listed in tablc[T] 
The measurements were started after 1500 equilibration trajectories and are based on 5000 
equilibrated trajectories. 



and multiple time scale integration described in detail in Ref. [8]. The trajectory 
length r was set to r = 1/2 in all our runs. Our estimates of the plaquette 
integrated autocorrelation time r mt (P) quoted in tableware in units of r = 1/2. 
Note that our estimates of the autocorrelation times of quantities such as amps 
or a/ps are found to be substantially smaller, typically by a factor of 5-10, than 
those reported in the table for the plaquette. 

As discussed above, maximal twist is realised in our simulations by tuning 
m to obtain a vanishing PCAC quark mass am PC Ac at the smallest value a/i min 
of the twisted mass parameter a/i. From table [2] one can see that this condition 
has been numerically realised with good accuracy, which in this context means 
^pc ac (/^min) //-^min < ^Kqcd within statistical errors (oAqcd ~ 0.1 in our case). 
Once this is achieved, the (weak) /t-dependence of mpcAc, which is visible in fig. 
1(a), is an 0(a) cutoff effect that merely modifies the 0(a 2 ) artifacts in physical 
observables, as already mentioned in section 2.1. 

In order to make maximum use of the gauge configurations, we evaluate con- 
nected meson correlators using a stochastic method to include all spatial sources. 
The method involves a stochastic source (Z(2)-noise in both real and imaginary 
part) for all colour and spatial indices at one Euclidean time slice. By solving 
for the propagator from this source for each of the 4 spin components, we can 
construct zero-momentum meson correlators from any bilinear at the source and 
sink. Four inversions of the Dirac matrix per Euclidean time slice value are nec- 
essary, since we chose to use only one stochastic sample per gauge configuration. 
This "one-end" method is similar to that pioneered in Ref. [66] and implemented 
in Ref. [HZj. We also employ a fuzzed source [68J of the extent of 6 lattice spac- 
ings to enable studying non-local meson creation and destruction operators. This 
allowed us to obtain very stable effective masses and to confirm the extraction of 
the pion ground state. 

In general, we save a gauge configuration every second trajectory and analyse 
meson correlators as described above from a selection of different Euclidean time 
slice sources. To reduce autocorrelations, we only use the same time slice source 
every 8-10 trajectories. Our primary statistical error was obtained with the so 
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Figure 1: (a) PCAC quark mass ampcAC as function of a/i and (b) Sommcr parameter (ro/a) 
as functions of (a/i) 2 . The solid line in subfigure (b) represents a linear fit in (a/i) 2 to the data. 



called r-method as described in Ref. [69] and cross-checked with a bootstrap 
analysis and a jack-knife analysis of blocked data. For a detailed description of 
our error analysis we refer to a forthcoming paper of our collaboration 



3.2 Force parameter rg 

In simulations of the quenched approximation of lattice QCD, the Sommer pa- 
rameter ro [TO] with a value of 0.5 fm, was widely used to set the lattice scale. 
While (r Q /a) is measurable to good accuracy in lattice QCD simulations it has 
the drawback that its value in physical units is not known very well. Therefore, 
it becomes necessary to determine the scale using other quantities which are ex- 
perimentally accessible with high precision, such as m n , f n , ttik, /k or uik*- In 
fact, in this paper we attempt to determine the lattice scale by fitting xPT based 
formulae to our precise data for / PS and m PS , using the physical values for m n 
and f n as inputs. From this analysis, we obtain a value of the lattice spacing 
which is 10% lower than the value obtained by setting r = 0.5 fm. 



Our results for (ro/a) are reported in tableland plotted in figure 1(b) Within 
the current errors the mass dependence of this quantity appears to be weak. Since 
r is a pure gauge quantity, it should be a function of (a/i) 2 and indeed, a linear 



fit in (a/i) 2 describes the data rather well as shown in figure 1(b) From the fit 
we obtain a value for r /a = 5.22(2) at the physical point, where a/i = a/j, n (see 
below), as also quoted in table [U 
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3.3 fps and mps as a function of the quark mass 

The charged pseudo scalar meson mass amps is as usual extracted from the time 
exponential decay of appropriate correlation functions [E In contrast to pure 
Wilson fermions, for maximally twisted mass fermions an exact lattice Ward 
identity allows to extract the (charged) pseudo scalar meson decay constant / PS 
from the relation 

/ps = 1(01^(0) |tt) I, (4) 

m PS 

with no need to compute any renormalisation constant since Zp = [2]. We 
give our results for m PS and / PS in table [21 

We now discuss whether the continuum x?T formulae can reproduce the data 
in table [2] for am P s and a/ P s- In our xPT based analysis, we take into account 
finite size corrections because on our lattices at the lowest and next-to-lowest 
/i-values they turn out to affect am P s and, in particular, a/ P s in a significant 
way. We have used continuum to describe consistently the dependence of 
the data both on the finite spatial size (L) and on /x. 

We fit the appropriate (Nf = 2) %PT formulae [HJ [72] 

2 

[l+eiog(2S Ax/A^] , (5) 

/ PS (L) — F [1 — 2Cgx (A)] [l - 2£ log(25 /// A 2 )] , (6) 
to our raw data for m P s and fps simultaneously. Here 

f = 2B li/(4ttF) 2 , A = y/2B ofJ ,L^ • (7) 

The finite size correction function gi(X) was first computed by Gasser and 
Leutwyler in Ref. [71] and is also discussed in Ref. [72] from which we take our 
notation (except that our normalisation of f n is 130.7 MeV). In Eqs. (jSj) and (151) 
NNLO ;\RT corrections are assumed to be negligible. The formulae above depend 
on four unknown parameters, Bq, F, A3 and A4, which will be determined by the 
fit. 

We determine afi n , the value of a fx at which the pion assumes its phys- 
ical mass, by requiring that the ratio [^/[mp s (L = oo)]// PS (L = 00)] takes 
the value (139.6/130.7) = 1.068. From the knowledge of afi n we can evaluate 
i 3)4 = log(A3 4 /m 2 ) and using f n the value of the lattice spacing a in fm. 

In order to estimate the statistical errors affecting our fit values we generate 
at each of the /i- values 1000 bootstrap samples for m P s and / P s extracted from 
the bare correlators, blocked with block-size of 32. For each sample (combining 
all masses) we then fit mp s and / P s simultaneously as a function of /i. From 
the 1000 fits we obtain 1000 bootstrap samples for 2aB Q , aF, log(a 2 Ag 4 ), 

2 Concerning results on the neutral pion mass, am-pg, see section (|3.4p . 



PS 



(L) = 2B li 
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and ^4, respectively, from which we compute the corresponding error estimates, 
taking in this way the statistical correlation between /pg and mps into account. 

For our lightest four /x- values, we find an excellent fit to our data on /pg and 
mps (see figures |2] and [3]) . The fitted values of the four parameters are 

2aB = 4.99(6) , 

aF = 0.0534(6), 

log(a 2 A 2 ) = -1.93(10), (8) 

log(a 2 A 2 ) = -1.06(4) . 



Our data are clearly sensitive to A 3 as visualised in figure 2(a). We obtain 



afi n = 0.00078(2), l 3 = 3.65(12), l 4 = 4.52(06) (9) 

which compares nicely with other determinations (for a review see Ref. [73]). 
Including also our results from a/i = 0.0150 in the fit gives an acceptable de- 



scription of mp s but misses the data for / PS , as shown in figures 2(b) and [3J 
Note, however, that in Eqs. (0 ED, and thus in the fit results (jSJ EJ), a number of 
systematic errors as discussed below are not included. 

The values presented here should hence be taken as a first estimate, the valid- 
ity of which will be checked in the future. Nevertheless, the statistical accuracy 
we are able to achieve implies that there is a very good prospect of obtaining 
accurate and reliable values for the low-energy constants from Wilson twisted 
mass fermion simulations. 

Since we have obtained an excellent description of our pseudo scalar mesons, 
we can use our fit to extract the lattice spacing. Based on the physical value of 
U, we get 

a = 0.087(1) fm. (10) 

Using the value of ro/a reported in table [T], this lattice calibration method yields 
r = 0.454(7) fm. 

We now discuss the possible sources of systematic error. Our analysis is based 
on lattice determinations of properties of pseudo scalar mesons with masses in 
the range 300 to 500 MeV on lattices with a spatial size slightly above 2 fm. 
Systematic errors can arise from several sources: 

(i) Finite lattice spacing effects. Preliminary results at a smaller value of the lat- 
tice spacing that were presented in Refs. [32], [33] suggest that 0(a) improvement 
is nicely at work and that residual 0(a 2 ) effects are small. 

(ii) Finite size effects. In order to check that next to leading order (continuum) 
xPT adequately describes these, we are presently performing a run at (3 — 3.9 
and ail = 0.004 on a 32 3 • 64 lattice. 

(iii) Mass difference of charged and neutral pseudo scalar meson. In the ap- 
propriate lattice xPT power-counting for our values of the lattice spacing and 
quark masses, i.e. a ~ /i ~ p 2 , one gets the order of magnitude relation 
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Figure 2: In (a) we show (amps) 2 /(a^) as a function of a^. In addition we plot the xPT fit 
with Eq. ^ to the data from the lowest four ^-values. In (b) we show (amps) 2 as a function 
of a[i. Here we present two %PT fits with Eq. ([5]), one taking all data points and one leaving 
out the point at the largest value a/i = 0.015. In both figures (a) and (b) we show finite size 
corrected (L — > oo) data points. 

(^ps) 2 _ ( m r>s) 2 = ^(^^-qcd) = 0(.P 4 )y fro m which it follows that to the order 
we have been working the effects of the pion mass splitting do not affect, in par- 
ticular, the finite size correction factors for m-ps and fp$. In spite of these formal 
remarks, it is possible, however, that the fact that the neutral pion is lighter 
than the charged one (by about 20% at a/i = 0.0040, see section (13.41) ) makes 
inadequate the continuum description of finite size effects adopted in the 
present analysis. This caveat represents a further motivation for simulations on 
larger lattices, which will eventually resolve the issue. 

(iv) Extrapolation to physical quark masses. We are assuming that at next 
to leading order for the Nf = 2 case is appropriate to describe the quark mass 
dependence of mp s and fps up to ~ 450-500 MeV. Our lattice data are consistent 
with this, but it would be useful to include higher order terms in the fits as 
well as more values of a/i to check this assumption. The effect of strange quarks 
in the sea should also be explored. 

3.4 Effects of Isospin Breaking 

In this section we report the results of some quantitative investigation of the 
effects of isospin breaking in the twisted mass formulation of lattice QCD at 
finite lattice spacing. This effect is expected to be largest in the mass splitting 
between the lightest charged and uncharged pseudo scalar mesons. A first analysis 
at a\i = 0.004, taking the disconnected contribution in the neutral channel fully 
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Figure 3: Wc show a/pg as a function of afi together with fits to xPT formula Eq. (JH|). 
We present two fits, one taking all data and one leaving out the point at the largest value 
a/i = 0.015. We show finite size corrected (L — > oo) data points. 

into account, shows that the uncharged pseudo scalar meson is about 20% lighter 
than the charged one. We obtain 



or, expressed differently, ^((mpg) 2 — (mp S ) 2 ) = c(a/r ) 2 with c = —4.5(1.8). 
This coefficient is a factor of 2 smaller than the value found in quenched investi- 
gations [2H]- Note that the uncharged pion being lighter than the charged one is 
compatible with predictions from lattice if the first order phase transition 

scenario is realised H7J HI] . For an investigation of isospin breaking effects in 
XPT see also Ref. p3]. 

The disconnected correlations needed for the 7r° meson are evaluated using a 
stochastic (Gaussian) volume source with 4 levels of hopping-parameter variance 
reduction [7S]. We use 24 stochastic sources per gauge configuration and evaluate 
the relevant propagators every 10-th trajectory. 

4 Summary 

In this letter we have presented results of simulations of lattice QCD with Nf = 2 
maximally twisted Wilson quarks at a fixed value of the lattice spacing a < 0.1 fm. 
We reached a pseudo scalar meson mass of about 300 MeV. The numerical 
stability and smoothness of the simulations allowed us to obtain precise results 
for the pseudo scalar mass and decay constant which in turn led to determine 
the low energy constants of the effective chiral Lagrangian. In particular, we find 
for the pseudo scalar decay constant in the chiral limit F = 121.3(7) MeV, and 
I3 = 3.65(12) and I4 = 4.52(6) where only statistical errors are given. 



am 



± _ 
PS — 



0.1359(7), 



am° PS = 0.111(11), 
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We do see effects of isospin breaking which are largest in the mass splitting of 
the neutral and charged pions and turn out to be about 20%. This is significantly 
smaller and opposite in sign than the corresponding splitting obtained in the 
quenched approximation. 

Tuning to maximal twist had to be performed on lattices of the same size as 
those used for the calculation of physical quantities. The reason for this is that 
we need to single out cleanly the one pion sector in order to impose the vanishing 
of the PCAC quark mass (Eq. (J3J)) without being affected by finite size effects 
or excited state contributions. Thus the tuning step itself is rather expensive. 
But it has to be done only once, as is the case for the determination of action 
improvement coefficients in other Wilson based approaches. Note, however, that 
with twisted mass fermions we do not have to compute further operator-specific 
improvement coefficients. 

The encouraging results presented here will be extended and checked by future 
simulations that will cover one coarser and one finer lattice spacing, double the 
statistics at one of our present simulation points {j3 = 3.9, a/i = 0.004) and go 
to a larger, 32 3 • 64, volume at the latter simulation point. In this way, we will 
be able to obtain results in the continuum limit, cross-check our autocorrelation 
times, improve our error estimates and control the finite size effects in order 
to check %PT predictions. The preliminary results presented in Refs. [32} [33] 
indicate a very good scaling behaviour already suggesting that automatic 0(a) 
improvement is indeed working well. 
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